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Abstract 

Whether in the Standard Model or beyond it, neutrinos contribute to the invisible 
decay mode of orthopositronium but practically not at all to that of parapositronium. 
Although this remark does not resolve the orthopositronium decay puzzle, it allows 
for upper bounds to be set on neutrino magnetic moments. 
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1. Positronium decay |2j has provided for a number of years a somewhat uncomfortable 
crack in this otherwise most beautiful edifice which is QED. Indeed, the two most precise 
experimental rates for the C-odd triplet orthopositronium (o-PS) decay[[| £3|, 

A^ xp = 7.0514 ± 0.0014 ^s" 1 , At P = 7.0482 ± 0.0016 /is" 1 , (1) 

are in flagrant conflict with the QED prediction |5], |6|, |2j, 

A? ED (PS -> 37) = 7.03831 ± 0.00005 /^s" 1 . (2) 

Note that the experimental results correspond, respectively, to a 9.4 a and a 6.2 o deviation 
from the theoretical expectation, and in relative terms, to the differences, 
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1.86 or 1.41) • 10" 3 . (3) 



In contradistinction, the C-even singlet parapositronium (p-PS) decay rate, of which the 
most precise measurement gives 0, 

A* xp = 7990.0 ± 1.7 fis' 1 , (4) 

compares better with the theoretical expectation [|]. [|. 

Af D (PS^2 7 ) = i^|^[l-"f5-^-2a 2 lna + ---l = 7989.5 /is" 1 . (5) 
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Recently however, a new measurement of the o-PS decay rate has been published^], 

A^ xp = 7.0398 ± 0.0025(stat.) ± 0.0015(syst.) /is" 1 , (6) 

in good agreement with the value in @, since it corresponds to a 0.5 a deviation from the 
theoretical prediction and a relative difference with it of 2.12-10 -4 . Clearly, an independent 
experimental confirmation of this beautiful result is desirable, before definitely concluding 
that the orthopositronium problem was indeed an experimental one. 

Over the years, many possible explanations for the discrepancy in the o-PS decay rate 
have been suggested and analysed, both theoretically and experimentally @. Among these, 
one may mention mirror photons ]10|], decay modes involving weakly coupled (pseudo) scalar 
bosons or even so-called invisible decay modes. However, none of these possibilities can 
accomodate the complete o-PS discrepancy above. For example, (pseudo)scalar decay 
modes typically lead to the following upper limit on the o-PS branching ratio[|], Pll - 

Br(o - PS -> 7 + X) < 10" 4 - 10~ 6 , (7) 

depending on the mass of the particle X, while the most stringent upper limit on the 
invisible branching ratio isfi~2l, |2|, 



Br(o - PS -> "invisible") < 2.8 ■ 10" 6 . (8) 

Clearly, given the remark in (|3]), these upper bounds exclude the possibility of explaining 
the o-PS problem in terms of such decay modes only. 
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Confronted with this difficulty, theorists have returned back[Q] to their calculations 
trying to establish the next correction^ in (a/n) 2 to the perturbative expansions || |], § 
used to determine the theoretical result in fl2[), 
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A? ED (PS - 3 7 ) = ^1 AZL^l [l + (-10.282±0.003)^ + ^ 2 lna+C (°) 

(9) 

However, in order to explain the observed discrepancy completely in terms of the next 
order calculation only, a coefficient C of the order of C ~ 250 ± 40 is required 0. Even 
though such a large coefficient is not to be excluded necessarily, such a value for the 
(a/n) 2 correction is difficult to contemplate. Incidentally, note that a value of C ~ 40 only 
is required by the recent new experimental result || given in (|6]). 

Curiously enough, although the above discrepancy seems to affect only o-PS decay, 
none of the suggested mechanisms has tried to exploit this possible hint towards an expla- 
nation. One such instance is that of invisible decay modes involving neutrino-antineutrino 
pairs. Indeed, even when slightly massive, (anti)neutrinos emitted in positronium decay 
would be essentially (right)left-handed, since this is their chirality in the Standard Model 
(SM) whether there exists physics beyond it or not. Therefore, momentum and angular 
momentum conservation implies that neutrino decay of p-PS is strongly suppressed — by a 
factor (m u /m e ) 2 — or vanishing altogether in the Standard Model, whereas even in the SM 
o-PS does decay via neutrino pair emission, albeit with the typical small rate of a weak 
interaction process. 

Obviously, in spite of the fact that this mechanism does indeed distinguish the decay 
of the two positronium hyperfine levels, it is expected to be far too small to explain the 
discrepancy for the o-PS rate, the relative weak interaction contribution having to be of 
the order (G ' pm^ 2 j ' o? ~ 2.4- 10~ 17 , Gp being Fermi's coupling constant. Nevertheless, the 
neutrino decay channel open essentially only to o-PS may be enhanced given additional 
couplings of the neutrinos. The possibility explicitly considered in this letter is that of 
non vanishing magnetic moments for Dirac neutrinos. Indeed, C-odd o-PS decay may then 
proceed through a single virtual photon which couples to the electron-positron pair and 
decays into a neutrino pair, while this additional decay channel is forbidden for the C-even 
p-PS state by charge conjugation invariance. Therefore, non vanishing neutrino magnetic 
moments would induce invisible decay modes of o-PS, but not of p-PS. In this way, given 
the measured o-PS decay rate, it is possible to establish upper bounds for neutrino magnetic 
moments. Such an analysis is the purpose of this letter^. 

Present experimental upper limits on neutrino magnetic moments are[| as follows! Id], 



H Ve < 1.08 ■ 10" 9 n B , /i. < 7.4 • 10" 10 imb , /V < 5.4 • 10~ 7 \i B , (10) 



fM B being the Bohr magneton, fi B — eh/2m e . Let us also recall the value [|l^, [16| of the 
magnetic moment of a massive Dirac neutrino in the SM, 

3eG F m e m u 19 (m v c 2 \ 

" y= 8W2 =3 - 2 " 10 \TeV r* ■ (U) 



1 It may be shown |]l4| |) that 57, 77, etc decay modes of o-PS cannot explain the discrepancy either. 

2 To the authors' knowledge, the only other study of PS neutrino decay appears in Ref. Ea]. However, 
that work addressed rather the radiative neutrino decay mode PS — > vVj only for massless neutrinos in 
the SM, and did not consider the possibility of non vanishing magnetic moments. 

3 Even though astrophysical or cosmological constraints lead to more stringent upper bounds |l6| |T7j ], 
such limits are model dependent and are not included here. 
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2. Let us first consider neutrino decay of positronium within the SM alone. For the purpose 
of applications beyond the SM however, neutrinos are taken to have a non vanishing Dirac 
mass already, but the possibility of flavour mixing will be ignored in this paper. The other 
assumptions entering the analysis are, on the one hand, that the couplings of these massive 
neutrinos to the W and Z are the usual ones as given by the SM, and on the other hand, 
that possible neutral Higgs contributions are not included since they are expected to be 
extremely small for light neutrinos. Relative to W and Z exchange diagrams, the neutral 
Higgs exchange amplitude is reduced by a factor m e m v /m\, a very small ratio indeed. 
Finally, two further approximations are effected; on the one hand, the Ps binding energy 
is neglected compared to the e~e + total rest- mass energy 2m e c 2 , namely the electron and 
positron are taken to annihilate at rest, and on the other hand, products of ratios of the 
electron or neutrino masses to the W and Z masses which appear in W and Z propagators 
are taken to be negligible as compared to unity. 

These assumptions having been stated, only two amplitudes may contribute to the 
neutrino decay process. On the one hand, there is the Z exchange diagram which con- 
tributes for all neutrino flavours]] vi (£ = e, //, r). On the other hand, there is the charged 
current amplitude which contributes only to the electron flavour neutrino channel u e via 
W exchange. A straightforward calculation then leads to the following contributions to 
the total decay rates. For the singlet state, one finds, 




e,fi,r 



(12) 
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where this last contribution follows from the interference of the W and Z exchange dia- 
grams. 

Similarly for the triplet state, one finds, 
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and for the interference contribution, 
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where 8\y is the usual weak mixing angle for neutral currents. 

4 For the sake of the analysis, three neutrino flavours whose mass is less than the electron mass are 
assumed to be involved in the decay process. Obviously, the argument would go through for an arbitrary 
number of neutrino flavours allowed by the decay kinematics. 
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A few remarks are in order. First, as was expected for reasons of the (right) left-handed 
chirality of (anti)neutrino couplings both to charged and to neutral currents, contributions 
to the singlet decay rate are all suppressed by a factor m 2 /m 2 . Second, even for the triplet 
state for which this chiral suppression is not effective, Z exchange contributions involve the 
factor (1 — 4 sin 2 9w = 0.0724) leading nevertheless to some suppression as well. And third, 
the relevant factor involving the weak coupling constant is (m 2 G 'p / 'ft a/2) 2 = 4.7-10 -25 , to be 
compared to the factors a 2 = 5.33-10 -5 and (2(7r 2 — 9)a 3 /97r) = 2.39-10 -8 relevant to the 27 
and 37 decays of p-PS and o-PS, respectively. Consequently, neutrino decay of positronium 
in the SM is very much suppressed, beyond the reach of any experiment at present. For the 
sake of the illustration, in the case of a massless neutrino [I 7^ e) one finds for example 
l/X { T Z \u e J^ e) ~ 5.14 -10 5 years, to be compared to l/Ay ED (37) ~ 142 ns! Incidentally, 
note that factors of the form (1 — am 2 /m 2 ) with (a ~ 1), do not differ significantly from 
unity for neutrino masses less than say, a fifth of the electron mass. 

3. Let us now extend the above analysis to include the possibility of neutrino ma- 
gnetic moments. As is well known [O, this requires massive neutrinos of the Dirac type, 
hence necessarily new physics beyond that of the Standard Model. The effective magnetic 
moment coupling of Dirac neutrinos is of the form, 



;is) 



where F af3 is the usual electromagnetic field strength F af3 = d a Ap — d@A a , and \i v the 



neutrino (anomalous) magnetic moment^]. Given this coupling, positronium decay then 
proceeds via one more amplitude in addition to those above, namely through the single 
photon annihilation channel e~e + — > 7 — > viPl. In the decay rate, this photon amplitude 
also interferes with the charged and neutral current amplitudes of the Standard Model. 

For the singlet state, one then finds the identically vanishing contributions, 



, e = 



e,/i,r 



(19) 

(20) 
(21) 



as is indeed required by charge conjugation invariance of electromagnetic couplings. On 
the other hand for the triplet state, one obtains, 
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(23) 



(24) 



In these expressions, \ib = eh/2m e is the Bohr magneton. Note that the (jZ) and (jW) 
interference contributions involve directly the ratio m u /m e , in contradistinction to the pure 
(77) contribution. The reason for this result is that the magnetic neutrino coupling implies 
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a spin flip whereas the W and Z (anti)neutrino couplings are purely (right)left-handed; 
angular momentum conservation thus requires one insertion of the neutrino mass vertex 
operator m u ip u ip u . Hence, given small neutrino masses as compared to the electron mass, 
it is essentially the (77) contribution which dominates over all other magnetic moment 
contributions, provided the factor a.\i v j\iB is of the order of the effective weak coupling 
{m L e G F l'K\pi) or larger. 

4. Given the above results, the total decay rates for both singlet and triplet states 
into neutrino flavours are obtained as, respectively, 
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and, 
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In particular, in the limit of vanishing neutrino masses, these expressions reduce to, 

A s (p - PS - i/P) = , (27) 
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where = 3 is the number of light neutrinos. Note that these two expressions confirm the 
announced result, namely the fact that in practice neutrino disintegration of positronium 
is a decay channel open essentially only to the triplet hyperfine state. 

The SM contribution from W and Z exchange relative to the QED 37 decay rate 
being 7.2 • 10~ 17 , let us first assume that the neutrino magnetic moment contributions 
to Xt (o — PS — > vv) are dominant, namely let us restrict the discussion now to the first 
term only in (p8| ). Compared to the QED decay rate A^ ED (PS — > 37) in @, the quantity 
relevant for the confrontation with experimental results is thus, 
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Taken at face value, if this ratio were to explain completely the discrepancy expressed 
in relative terms in one would require the following limit, 



J2 »\ ^3-10-Vb 



(30) 



-e,fj,,T 



On the other hand, if the corresponding limit stemming from the more recent measurement 
in (pD is used in the same manner, one derives, 



E 



(31) 
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However, since the upper limit on the invisible branching ratio in (H) already excludes such 
possibilities, it is more sensible to use the constraint on this mode to set the upper bound, 



J2 l4 t <1.2-KT 4 /ib . (32) 



=e,^t,r 



Clearly, these numbers are not competitive with the experimental limits on neutrino ma- 
gnetic moments in fllTf). Even for the least stringent upper bound which applies in the case 
of the neutrino v T , \l Vt < 5.4 • 10~ 7 fi B , the corresponding branching ratio is already, 

Br(o - PS -> u T v;) < 5.5 • 10 -11 , (33) 

namely much less than present experimental upper limits on branching ratios for invisible 
decay modes or for the relative experimental deviations from the QED prediction. 

As a matter of fact, it is instructive to consider the situation when the magnetic 
moment contribution to the neutrino branching ratio is comparable to that of the charged 
and neutral currents, namely when a^/f^s = m^Cp /7r\/2- Indeed, the corresponding 
magnetic moment value, 

\i v = 9.39 ■ 10 -11 fi B , (34) 

is of the same order of magnitude as the best experimental upper limit in ([H]) established 
for the muon neutrino. A posteriori, this is not surprising, since these experimental limits 
are determined from neutrino scattering experiments designed to be sensitive to processes 
whose strength is typical of weak interactions. Hence, present limits on neutrino magnetic 
moments necessarily correspond to values which render their contributions comparable to 
those of the ordinary charged and neutral electroweak currents. In the present case, the 
magnetic moment value in (|34]) for a single neutrino contributes the invisible branching 
ratio the quantity 1.7 • 10~ 18 , which is indeed close to the neutrino branching ratio in the 
SM of 7.2 • 10~ 18 , namely the contribution of the second term in (p8|). 

In fact, using the upper limits in (|10|) as values for the neutrino magnetic moments 
together with the result in (|28f) , one obtains the following total neutrino branching ratio, 

Br(o - PS -> vv) = 5.6 • 10 -11 , (35) 

which is thus dominated by the magnetic moment contribution of the v T neutrino. 

However, the value in fl3~5] ) is larger than upper bounds recently established |]ST| for 
the branching ratio of "exotic" positronium decays on the basis of primordial nucleosyn- 
thesis. In fact, since the electroweak charged and neutral current contributions relevant to 



the thermal equilibrium of three light neutrinos are already included [22] in the standard 



cosmological model, the "exotic" contributions to be considered here are solely those stem- 



ming from the neutrino magnetic moments. Moreover, the associated couplings (|18|) being 
dimension five operators, it is actually the upper bound of 2 • 10~ 17 on the branching ratio, 



associated to case B) of Ref.pi|], which is relevant in our case and which thus applies only 
to the total contribution of neutrino magnetic moments. Correspondingly, this upper limit 
implies, 

J E < S3.2.1(T 1( Vb • (36) 

Y £=e,fi,T 

Incidentally, note that this upper bound value is also typical of astrophysical constraints 
on neutrino magnetic moment s [|lq] . 
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5. In conclusion, even though neutrino decay of positronium in the Standard Model 
and beyond it does indeed contribute to invisible decay modes of the triplet orthopositro- 
nium state, but practically not to the decay of the singlet parapositronium state, present 
experimental limits on neutrino magnetic moments imply that such processes cannot pro- 
vide even for a partial contribution towards a resolution of the possible orthopositronium 
lifetime puzzle [Q], which however, may have been resolved by a recent new measurement . 
Moreover, these limits also establish that the observation of such decay modes is far be- 
yond present experimental capabilities. Nevertheless, when considered together with recent 
argument s [pT|l based on primordial nucleosynthesis, the orthopositronium neutrino decay 
mode implies upper bounds on neutrino magnetic moments which are more stringent than 
present experimental limits. 

Although the analysis assumes massive Dirac neutrinos without flavour mixing, it 
should be clear that similar conclusions would apply more generally for transition magnetic 
moments as well, whether the massive neutrinos are Dirac or Majorana spinors with flavour 
mixing. 

Note also that an analogous discussion could be developed were neutrinos to carry 
electric charges |2"3| eQ v . Under such circumstances, the associated couplings would imply 
contributions to positronium decay such as those above in which the ratio cl\i v Ihb is 
replaced by the factor aQ u , thereby leading to similar conclusions. For example in the 
case of the triplet decay rateQ, the modulus squared photon amplitude leads to, 



A^ 7 (ye; Q 



1 3 



{aQv 
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1 m 
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+ 2 



mi 



(37) 



Therefore, the same types of upper bounds as those derived above for the ratios fi u /fiB 
would apply to the neutrino electric charges Q v . However, these constraints are not com- 
petitive with existing limitspS. In this respect, it is of interest to remark that based on 



the anisotropy of the microwave background, an upper bound of \Q V \ < 4.8 ■ 10 was 



recently obtained |24|], which is thus many orders of magnitude more stringent than any 
existing limit on neutrino magnetic moments. 



It is a pleasure to thank Prof. Jean Pestieau for his interest in this work and for 
discussions. 



5 Here again, contributions to the singlet rate vanish by charge conjugation invariance. 



7 



References 

[1] For reviews on positronium, see for example, 
A. Rich, Rev. Mod. Phys. 53 (1981) 127; 

A. P. Mills, Jr. and S. Chu, in Quantum Electrodynamics, ed. T. Kinoshita (World 
Scientific, Singapore, 1990), pp. 774-821. 

[2] The present letter shall not attempt a review of the o-PS lifetime problem, covered 
for example in, 

M.I. Dobroliubov, S.N. Gninenko, A.Yu. Ignatiev and V.A. Matveev, Int. J. Mod. 
Phys. A8 (1993) 2859; 

V.V. Dvoeglazov, R.N. Faustov and Yu.N. Tyukhtyaev, Mod. Phys. Lett. A8 (1993) 
3263. 

[3] C.I. Westbrook, D.W. Gidley, R.S. Conti and A. Rich, Phys. Rev. Lett. 58 (1987) 
1328; (E) ibid 58 (1987) 2153; Phys. Rev. A40 (1989) 5489. 

[4] J.S. Nico, D.W. Gidley, A. Rich and PW. Zitzewitz, Phys. Rev. Lett. 65 (1990) 1344. 

[5] A. Ore and J.L. Powell, Phys. Rev. 75 (1949) 1696; 

W.E. Caswell, G.P Lepage and J.R. Sapirstein, Phys. Rev. Lett. 38 (1977) 488; 

W.E. Caswell and G.P. Lepage, Phys. Rev. A20 (1979) 36; 

M.A. Stroscio, Phys. Rev. Lett. 48 (1982) 571; 

G.S. Adkins, Ann. Phys. (N.Y.) 146 (1983) 78; 

G.S. Adkins, Phys. Rev. A27 (1983) 530; ibid A31 (1985) 1250; 

G.S. Adkins, A.A. Salahuddin and K.E. Schalm, Phys. Rev. A45 (1992) 3333, 7774. 

[6] LB. Khriplovich and A.S. Yelkhovsky, Phys. Lett. B246 (1990) 520. 

[7] A.H. Al-Ramadhan and D.W. Gidley, Phys. Rev. Lett. 72 (1994) 1632. 

[8] P.A.M. Dirac, Proc. Cambridge Phil. Soc. 26 (1930) 361; 
I. Harris and L. Brown, Phys. Rev. 105 (1957) 1656; 
W.E. Caswell and G.P. Lepage, Phys. Rev. A20 (1979) 36; 
Y. Tomozawa, Ann. Phys. (N.Y.) 128 (1980) 463. 

[9] S. Asai, S. Orito and N. Shinohara, Phys. Lett. B357 (1995) 475. 

[10] S.L. Glashow, Phys. Lett. B167 (1986) 35; 

E.D. Carlson and S.L. Glashow, Phys. Lett. B193 (1987) 168; 
S.N. Gninenko, Phys. Lett. B326 (1994) 317. 

[11] S. Asai, K. Shigekumi, T. Samuki and S. Orito, Phys. Lett. B323 (1994) 90; 

T. Maeno, M. Fujikawa, J. Kataoka, Y. Nishihara, S. Orito, K. Shigekuni and 
Y. Watanabe, Phys. Lett. B351 (1995) 574; 

T. Mitsui, K. Maki, S. Asai, Y. Ishisaki, R. Fujimoto, N. Muramoto, T. Sato, Y. Ueda, 
Y. Yamazaki and S. Orito, Eur. Phys. Lett. 32 (1996) 111. 

[12] T. Mitsui, R. Fujimoto, Y. Ishisaki, Y. Ueda, Y. Yamazaki, S. Asai and S. Orito, 
Phys. Rev. Lett. 70 (1993) 2265. 



8 



[13] See for example, 

P. Labelle, G.P. Lepage and U. Magnea, Phys. Rev. Lett. 72 (1994) 2006; 
M. Samuel and G. Li, Phys. Lett. B331 (1994) 114; 

A.P. Burichenko, Phys. At. Nucl. 56 (1993) 640 (Yad. Fiz. 56 (1993) 123); 

LB. Khriplovich and A.I. Milstein, J. Exp. Theor. Phys. 79 (1994) 379; 

A.P. Burichenko and D.Yu. Ivanov, Vacuum polarization corrections to the positronium 

decay rate, Budker preprint BINP 94-64 (November 1994), [hep-ph/9411322 . 

[14] G.P. Lepage, PB. Mackenzie, K.H. Streng and P.M. Zerwas, Phys. Rev. A28 (1983) 
3090; 

G.S. Adkins and F.R. Brown, Phys. Rev. A28 (1983) 1164. 
[15] W. Bernreuther and O. Nachtmann, Z. Phys. Cll (1981) 235. 
[16] Particle Data Group, Phys. Rev. D50 (1994) 1193. 

[17] J. Bernabeu, S.M. Bilenky, F.J. Botella and J. Segura, Nucl. Phys. B426 (1994) 434. 

[18] B.W. Lee and R.E. Shrock, Phys. Rev. D16 (1977) 1444; 

K. Fujikawa and R.E. Shrock, Phys. Rev. Lett. 45 (1980) 963. 

[19] See for example, 

R.N. Mohapatra and P.B. Pal, Massive Neutrinos in Physics and Astrophysics (World 
Scientific, Singapore, 1991). 

[20] C. Itzykson and J.-B. Zuber, Quantum Field Theory (McGraw-Hill, New York, 1980), 
pp. 66-67. 

[21] R. Escribano, E. Masso and R. Toldra, Phys. Lett. B356 (1995) 313. 
[22] See for example, 

S. Weinberg, Gravitation and Cosmology (Wiley & Sons, New York, 1972), pp. 469- 
609. 

[23] For a recent discussion and references, see for example, 

AYu. Ignatiev and G.C. Joshi, Phys. Rev. D51 (1995) 2411. 

[24] S. Sengupta and P.B. Pal, Phys. Lett. B365 (1996) 175. 



9 



